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The third post-Newtonian approximation to the general relativistic dynamics of two point-mass 
systems has been recently derived by two independent groups, using different approaches, and 
different coordinate systems. By explicitly exhibiting the map between the variables used in the two 
approaches we prove their physical equivalence. Our map allows one to transfer all the known results 
of the Arnowitt-Deser-Misner (ADM) approach to the harmonic-coordinates one: in particular, it 
gives the value of the harmonic-coordinates Lagrangian, and the expression of the ten conserved 
^T>' quantities associated to global Poincare invariance. 
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I. MOTIVATION 



m 

^ ■ Binary systems made of compact objects (neutron stars or black holes) are the most promising sources for the 
upcoming ground-based network of interferometric gravitational wave detectors LIGO/VIRGO/GEO. Because of 
their higher signal-to-noise ratio, the first detections are likely to involve massive binary black-hole systems, with 
total mass mi + m<i > 30 Mq. Such systems emit most of their useful signal at the end of their inspiral phase, near 
the last stable (circular) orbit. This makes it very important to have the best possible analytical control of the general 
' relativistic dynamics of two-body systems. 

For many years the equations of motion of binary systems have been known only up to the 5/2 post-Newtonian 
(2.5PN) approximation 0-0]. Recently, Jaranowski and Schafer j|J^] and Damour, Jaranowski, and Schafer jl0 11 



q-i succeeded in deriving the third post-Newtonian (3PN) dynamics of binary point-mass systems within the canonical 
formalism of Arnowitt, Deser, and Misner (ADM). More recently, Blanchet and Faye JlJ,[l3| succeeded in deriving the 
^Jy 3PN equations of motion of binary point-mass systems in harmonic coordinates relying on an independent framework. 
• • . The purpose of this paper is to compare and relate these two sets of results. 

For the present investigation it is sufficient to consider the conservative part of the dynamics, i.e. we shall drop 
the dissipative 2.5PN part which is the leading order radiation-damping level. This has the advantage that for the 
remaining part there exists an autonomous Hamiltonian (depending only on particle variables) and a conserved energy. 
The non-autonomous parts of the ADM Hamiltonian, up to the 3.5PN level (which is the next order radiation-damping 
level after 2.5PN), are given in Ref. Q. 



II. REGULARIZATION AMBIGUITIES 



Before tackling the comparison between the two sets of 3PN results several remarks are in order. First, let us 
emphasize that both approaches to the 3PN dynamics have found that the use of Dirac-delta-function sources to model 
the two-body system causes the appearance of both badly divergent integrals and badly defined "contact terms" , which 
(contrary to what happened at the 2.5PN J|,|| and 3.5PN |lj] levels) cannot be unambiguously regularized. More 
precisely, when Refs. |8||^] derived the relative- motion 3PN ADM Hamiltonian -ff (x, p), in the center-of-mass frame of 
the binary, they introduced two arbitrary dimensionless parameters, u>k{= ^kinetic) and uj s (= w s tatic), to formalize the 
presence of irreducible ambiguities in the regularization of the Hamiltonian. The regularization ambiguity parameter 
ujk concerned a momentum-dependent contribution cx G 3 c~ 6 (p 2 — 3(n • p) 2 )r~ 3 , while to s concerned a momentum- 
independent contribution oc G 4 c~ 6 r~ 4 . Ref. Jll]], on the one hand generalized the work of || by deriving the 3PN 
ADM Hamiltonian H(xi, xa, Pi, P2) in an arbitrary reference frame, and, on the other hand, proved that u>k was 
uniquely determined to have the value u>k = 41/24 by requiring the global Poincare invariance of the 3PN dynamics 



1 



(see Ref. [ fl5| for details of why uik is not fixable in the center-of-mass frame). Therefore, finally, the 3PN ADM 
Hamiltonian|j contains only one regularization ambiguity (parametrized by the dimensionless parameter u> s ): 



H(-K a ,p a ;uj s 



m a c 2 + iJ N (x a ,p a ) 



1 



#1PN(X , p a ) + ^-i?2PN(x Q , Pa) 



^3PN ( X oj Pa) 



G 4 m^m|(mi + 7712) 
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(2.1) 



See Ref. |Tl[] for the explicit expressions of iJ^p^ , and of the well known lower-order (Newtonian, 1PN, and 2PN) 
contributions. 

On the other hand, Ref. fl3|| , in deriving the 3PN equations of motion in harmonic coordinates, introduced four 
arbitrary dimensionfull parameters s\, S2, r[, and r 2 (with dimensions of length; s±, S2 correspond to the intermediate 
ADM regularization length scales h, h discussed below, which disappear in the final i?3PN). In addition, though these 
authors developed some formal generalization of the theory of distributions to deal with the badly divergent integrals 
appearing at 3PN |l6| ], they could not prove the uniqueness of their prescriptions, and, in fact, they used two different 
prescriptions, the most recent of them [ p"3"[ introducing a new dimensionless arbitrary parameter K . (However, they did 
prove that both versions of their regularization prescriptions finally lead to gauge-equivalent equations of motion.) In 
jl7j these authors also introduced a modification of their regularization procedures, aimed at yielding ("in principle" ) 
Lorentz-invariant equations of motion. The net result of using the set of regularization recipes developed in |l6|,|l7|,[l3| 
is the derivation of 3PN two-body equations of motion, in harmonic coordinates, which depend on the five parameters 
s ii s 2, r'2, and K and which are, generically, neither Lorentz- invariant nor deducible from an action (because 
they do not lead to a conserved energy as any autonomous action-based equations of motion would). Then, the 
authors of |l3[ ] impose the triple requirement of: (i) Lorentz invariance, (ii) existence of a conserved energy, and (iii) 
polynomiality in m\ and 7712. They show that: (i) uniquely determines K to have the value K = 41/160, and (ii) 
imposes one constraint relating the four length scales s\, S2, r[, and r 2 , namely 
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(2.2) 



Note that, when they use their older version of their regularization prescriptions, the rational number appearing in 
Eq. (2.2) becomes —159/308. By further imposing the requirement (iii), they conclude that the two length scales 
Si, S2 can be expressed in terms of the two other scales r[, r' 2 , and of a new dimensionless parameter A, through 
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(2.3a) 
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(2.3b) 



Finally, the 3PN equations of motion for the harmonic coordinates y a {t) of the two point masses contain three 
regularization ambiguities (parametrized by the two scales r[, r' 2 and the dimensionless parameter A) and have the 
form y a = A a (y D ,v ), where v a = y a , with 

Aa(y b ,v b ) = A oN (y 6 ,v 6 ) + -^A alPN (y b ,v b ) + -i-A a2 PN(y , v b ) + -^A a 2 . 5 PN(y&, v b ) 



Ai°3 ) PN (y b ,v ) +ln A« PN (rh 2 ,vi 2 ) + In A« PN (r? 2 , v 12 ) + A A^ PN (r? 2 ) 



(2.4) 



1 Note that we are considering here the ordinary 3PN Hamiltonian, obtained (following a result of ju]]) by a well-defined shift 
of phase-space coordinates, designed to reduce the higher-order Hamiltonian _ff3PN(x a , p a , x a , p a ) defined by eliminating the 
field variables hjj 1 , hjj 1 in the "Routh functional" J?3PN(x a , p a , hjj 1 , hf^) of ||. 

2 At least if one follows in using the new "correct" derivative involving the parameter K. 
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where r\ 2 = yi — y2 and V12 = = v i — v 2 denote the (harmonic) relative position and velocity, respectively. It 
was, however, shown in [H| that the ambiguities linked to r[, r' 2 can be gauged away, so that the physical ambiguity of 
the harmonic equations of motions is described by the sole (dimensionless) parameter A. For simplicity, we shall work 
here with the equations of motion explicitly displayed in |13| w hich , in fact, corresponds to their older regularization 
prescription [with 783/3080 being replaced by —159/308 in Eq. (|2.2|)]. See Eq. (7.16) of [jl3| for the explicit expression 
of the 3PN contributions to the harmonic equations of motion (as well as of the well known lower-order contributions 

|0,||). We shall only note here the fact that A^gp N and A^^p^ depend only on the relative positions and velocities, 
and that the A-term reads (for a = 1; n 1 ^ = r^/r^) 

AA (3) f M _ 44 GVml(mi+m 2 ) h 

AA 13PNV r 12j - ~y A ^3Tj5 n !2- ^• t> ) 

Even before any detailed calculation, it is clear that this A-contribution derives from a potential energy 

\T/( 3 )= n , G4m i w 2( m i + m 2) 

so that, if the two different 3PN dynamics can be shown to be somehow equivalent, the "harmonic" regularization 
ambiguity A must be related to the "ADM" one u s by — ^-A = uj s + const. 



III. ORIGIN OF REGULARIZATION AMBIGUITIES 

As the presence of the regularization ambiguities at the 3PN level is very striking^ and physically momentous, let 
us discuss in more detail the origin of the ambiguities present in the two approaches, and the differences between 
them. 

In the ADM approach, one computes a (spatially) global scalar quantity, the Hamiltonian _ff(x a ,p a ) of the system. 
Essentially^ the global scalar £f(x ,p a ) can be explicitly expressed as an integral over space of an integrand of the 
generic form 

H(x;x a ,p Q ) = ^ D) (x;x„,p a ) + W^x; x Q , p a ) + <9;L> 4 (x; x Q , p a ). (3.1) 

Here, Ti c is made only of "contact terms", i.e. of terms proportional to the delta-functions modelling the sources, say 
Tii D ^ — J2 a Safe, Xb)<5(x — x a ), where S a is constructed from field quantities, Ti^ is a "field-like" term, i.e. an "energy 
density" constructed from field quantities and distributed all over space, and the last term is a pure divergence, which 
formally gives a vanishing contribution^] to the integrated Hamiltonian. As indicated by the superscript notation, the 
explicit values of the "contact" and "field" terms depend on the choice of the divergence term diD 1 . In other words, 
we can, by "operating by parts", shuffle terms between Tic and TiP , at the price of changing D l . Note that, when 
so shuffling terms, one freely uses Einstein's field equations (with delta-function sources) and one assumes the validity 
of the usual rules of functional calculus,^] such as Leibniz' rule [di(AB) = (diA)B + A(diB)], and the commutativity 
of repeated derivatives (didjA — djdiA). 



3 Though it was anticipated in M, see pp. 107 and 116 there. 

4 After applying the double "reduction" process of eliminating the field variables and reducing the order of the Hamiltonian 
i@. 

It has been checked in the ADM approach that the "surface term at infinity" associated to diD 1 is not causing any ambiguity. 
Indeed, most pieces in <f dSiD 1 decay like some inverse power of r at infinity, while the ones which might be problematic (like 
the one associated to the 0(r) part of hf^j) have been explicitly shown to give a vanishing contribution to <f dSiD 1 . The 
ambiguities come only from the singular behaviour of the integrand near each particle, i.e. as the field point x tends to either 
xi or X2. 

6 In the explicit computations of the Hamiltonian done in Refs. pcj , and jll| one has always chosen D l 's such that all the 
terms in Ti^P^ contain only one derivative (or its equivalent) acting on the elementary fields (4>(2), 4>(A) > • • • ), so that there 
is no need to worry about using an improved distributional derivative. The distributional rule of differentiation of homogeneous 
functions described in Appendix B of || is, in fact, used only when gauging the ambiguities by computing the regularized value 
of J d 3 x(diD 1 ), as explained in the Appendix A of Jlo| , 
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The ambiguities in the determination of the value of H = J d 3 x7i come from two separate (but related) facts. 
First, the "contact" contribution 



ee J d 3 xH^ = J d 3 x^5 a (x,x 6 )5( 5 



is formally infinite because the (field-constructed) quantity 5 tt (x, X;,) is generically singular To give a 

meaning to H { C D) one must choose a specific regularization prescription to define the limit lim x ^ Xa S„(x, x&). Second, 
the "field" contribution = J d 3 xHj is also formally infinite because the integrand T~tf is generically too 

singular asx^x„ to be locally integrable. To give a meaning to H ^ one must choose a specific regularization pre- 
scription for such singular integrals. Finally for each choice of D l , one defines the regularized value of the Hamiltonian 

In Refs. [glU the following specific regularization prescriptions were adopted: (i) for contact terms 
[lim x ^ Xa 5 a (x, Xb)] rcg is defined as Hadamard's "partie finie" of 5 a (x, x&), Pf a 5 a (x, X&), defined in Appendix B of || as 
the angle-averaged finite term in the Laurent expansion of iS a (x a +r Q , x?,) in powers ofr a = \r a \ = |x — x a | (as r a — > 0), 
and (ii) for field terms the regularized value 7 rog of a singular integral I = j d 3 xF(x, xi, X2) is defined as follows. 
First, one regularizes separately the divergences near each particle, i.e. the integrals I a = J v d 3 xF(x. 7 xi, X2) where V a 
is a volume which contains x but not x&, with b ^ a. [Evidently, one can always decompose I = I\ + 1% + I C ompi with 
two local volumes Vi, V2 and a regular complement.] Second, each local integral, say I\ near particle 1, is regularized 
"a la Riesz" , i.e. by analytic continuation (AC) in t\ of 



Ji(ei) = ^ d3x (ji) F ( x ' x i> x 2) 



where n = x — xi| and where l\ is a certain length scale.[| Most integrands F lead to functions of e±, h(ei), which 
are analytically continuable into the complex ei-plane down to e\ = 0. In such a case this continuation AC ei ^o^i(ei) 
uniquely defines the regularized value of 1% . However, a limited subclass of "dangerous" integrals gives rise to a (simple) 
pole as e\ — > 0: ii(ei) = Zxfa 1 +\n(Ri/h)) + A u where Ri is an "infra-red" length scale associated to the choice of 
the local volume V\. For such integrals, one is naturally led (following the usual "minimal subtraction" prescription 
of quantum field theory) to defining the regularized value of I\(ei) as the limit of /i(ei) — Z\/ei as e\ — > 0, i.e. as 
J[ eg = Z\\vl{R\H\) + A\. Note that this regularization prescription has introduced one arbitrary length scale: the 
regularization length lx- [The -related infra-red length Ri is easily seen to cancel out in I = I\ + h + icompl-] 
However, as emphasized in Sec. IV of ||, a remarkable thing occurs in the explicit calculations of the 3PN ADM 
Hamiltonian: the combination of dangerous integrals appearing in i?3pjv is such that all pole terms exactly cancell: 
^2 Zi =0. In fact, one of the characteristics of the calculation of i?3PN in the ADM formalism is that one finds it 
much safer (and simpler) to regularize, at once, the full integral, rather than to try (as in the harmonic-coordinate 
calculation, |l3[| ) to give a separate regularized value for each individual contribution to the equations of motion. 
The global cancellation of the poles shows that the combination of dangerous integrals appearing in i?3PN is of a 
less dangerous type. A nice aspect of this cancellation is that, in the ADM approach, the two regularizing length 
scales li, I2 completely cancell and do not appear in the regularized final i?3PN- This does not mean, however, that 
the final result is unambiguous. Indeed, it was emphasized in ||,|l(| that the regularized value of i?3PN depends on 
the reshuffling of terms used to separate Ti in TL^ + T~t^ + d{D\ In other words, when operating by parts (which 
changes D\ and TC{ D ^ and Tif) the regularized value of = He + is found to change. In addition 

to this D'-dependent ambiguity, there is also the problem of the sensitivity of the contact contribution Hc°^ to the 
choice of prescription for defining the 'partie finie' of if? a (x,X6). It was emphasized in that the definition of the 

'Hadamard partie finie' Pf a becomes ambiguous at 3PN because it cannot be 'threaded' through a product of field 
functions, i.e. that, in general, Pf a (/i/2 ■ ■ •) ^ (Pfa/iXPfa/2) • • ■ [The prime, irreducible example of this ambiguity at 
3PN comes from the fact that Pf a (0^) ^ [Pf a (</>^)] 2 = [Pf a ((/>( 2 ))] 4 , where (f)(2) is the Newtonian potential.] 



7 The "Riesz" prescription explained in the Appendix of B of ^ looks different from what we explain here (because it does 
not separate the integration volume in Vi, V2 and the rest), but, as emphasized in it is equivalent to the logically simpler 
prescription that we summarize here. 

8 In actual calculations (see especially Appendix A of jl^) one monitors the changes in H^ D ^ = _ffi°' + by computing 
the term-by-term regularized value of the full algebraic expansion of the divergence term J d 3 x(diD l ). 
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The attitude of Refs. §-[k]] vis a vis these regularization ambiguities has been the following: (i) one must acknowl- 
edge their existence, because there exists, as yet, no convincingly unique extension of distribution theory allowing one 
to select a preferred regularized value, and (ii) however, one can analyze in detail the structure of these ambiguities 
and show that they can be parametrized by only two (dimcnsionlcss) parameters: lo^ and uj s . Indeed, after the pio- 
neering work which introduced these regularization ambiguity parameters, a systematic study of the ambiguities 
has been conducted in the Appendix A of [|10| (by exploring all the possible operations by parts, as well as the effect 
of having Pf(/i/2) ^ Pf(/i) Pf(/2))- This study confirmed the existence of only two regularization ambiguities As 
the most recent work in the ADM formalism has shown that the 'kinetic' ambiguity u>k was uniquely fixed by 
imposing global Poincare invariance, the final conclusion is, as indicated in Eq. fl2.l| ) above, that the ADM formalism 
introduces only one regularization ambiguity parameter: the 'static' ambiguity w s . 

It would take us too long to explain in detail why the harmonic-coordinate approach introduces more ambiguity 
parameters [four, (sj., s%, r[, r' 2 ), or five, (si, S2> f[, r' 2 , K), depending on the regularization prescription, instead of 
two, (ujk,0J s )] than the ADM one (see p6| , p^7| , ^3[ ) . Let us only make a short list of the most significant differences 
between the two approaches: (i) Blanchet and Faye regularize separately many independent singular contributions 
to the spatial derivative of the gravitational field instead of working with the full scalar Hamiltonian as a block, (ii) 
when computing their "elementary integrals" by analytic continuation they can (after contracting free indices) use 
the ordinary Riesz formula (instead of the generalized Riesz formula of Q , necessary to deal with the denominators 
oc (ri + r-2 + ri2) 7 that appear in the ADM Hamiltonian), (iii) they directly work with the full hierarchy of PN fields 
up to <?oo = — 1 + ■ • • + 2j7§/c 8 , while the ADM approach needs to work only with the contribution ifi^/c 6 to the 
"scalar" potential, (iv) they get two (gauge-related) ambiguities of "logarithmic" type (involving two arbitrary length 
scales) , and (v) they use a different coordinate system. It would be interesting to study whether a reworking of the 
harmonic-coordinates work along the more 'global', and more 'PN order reduced' lines of the ADM approach would 
not simplify their results and get rid of several of their ambiguities. 

IV. MATCHING THE TWO 3PN DYNAMICS 

We shall now show in detail that the two 3PN dynamics are equivalent modulo a suitable shift of particle vari- 
ables. Some time ago, Damour and Schafer || studied the link, at the 2PN level, between the ADM dynamics and 
the harmonic-coordinates (or DeDonder-coordinates) one. They explicitly constructed the map between these two 
descriptions of the dynamics. Let us emphasize that this map acts on the "motions", i.e. on the particle positions 
(and momenta or velocities) as functions of time. In other words, it gives either the transformation (with v a = y a ) 

ya(*)=Y„(x 6 (t),p a (t)) 1 (4.1a) 

v o (t)=V o (x6(0>P&(*)) s (4-lb) 
from the ADM variables (x(,,p(,) to the harmonic ones (y a ,v a ), or the inverse transformation 

x a (i) =X a (y 6 (t),v 6 (t)), (4.2a) 

P«(*)=P„(y6(t),v 6 (t)). (4.2b) 



All variables in Eqs. (4.1), (4.2) are taken at the same value for their (respective) time argument. As explaine d in 
it is always possible to express the looked for map in this form. One has to beware that the transformations (4.1) or 



(4.2) are not the direct restriction of a coordinate transformation, x ,fi = x^ + ^(x , [xi], [X2]) (where the brackets 
indicate functional dependence), to a field point x M on a particle world line, but that one must take into account the 
time-shift £° to transform the coordinate shift into the "motion" shift x'(i) — x(t ) [se e Eqs. (3) of ||] 



Among the two forms (4.1) or ([4.2]) we found that it is simplest to work with (4.1). Indeed, the necessary and 



sufficient conditions for (4.1) to map the ADM dynamics onto the harmonic one is easily seen to be simply 

{Y a ,H} = V ai (4.3a) 

{V a ,H} = {{Y a ,H},H} = A a (Y b ,V b ). (4.3b) 



9 Ref. |[o| made an attempt at lessening the sources of ambiguity by choosing a D % such that the contact terms Hi are 
absent. However, even in this 'preferred' presentation, gave rise to the two usual ADM ambiguities. 
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All functions entering Eqs. (4.3a) and ( 4.3b| ) (Y a ,H, V a ) are functions of the ADM phase space coordinates (xb,ph) 
The notation {•, •} denotes the usual Poisson bracket 



{A(x a , Pa ), B(x a , Pa )} =J2Y1 



dA dB OA dB 



dx a dpai dp ai dx a 



Fina lly H denotes the full 3PN Hamiltonian (2.1) while A a (yb,Vf,) denotes the harmonic equations of motion, Eq. 



(2^4). Note that Eq. ( 4.3a ) explicitly determines V a (xb, p b) in terms of Y a (xf,,pf,). Therefore, the problem of the 
mapping between H and A a is reduced to solving Eq. (4.3b) as an equation for the two unknown phase-space vectorial 
functions Y^x^p;,) and Y2(x;,,pb). We tackled this problem by the method of undetermined coefficients, i.e. by 
writing the most general expression for the PN expansion of Y a (xb,p&). We know that Y a differs from x a only at 
2PN order, i.e. 



Y a (x b , Pb ) = x a + lyf N (x b , Pb ) + ivf N (x b ,p fc ). 



(4.4) 



The explicit expression of Y^ PN was given in Ref. || (we write it here for a = 1; the expression for a = 2 being 
obtained by a simple relabeling 1 <-> 2): 



Yf N (x Q , Pa ) =Gm 2 < 



5 Pi 1 (ni2 • p 2 ) 



8m| 



Gm 1 (7_ 



1 ma 
4 mi 



ni2 



1 (ni2 ■ pa) pi _ 7 (ni2 ■ p 2 ) P2 

2 ma mi 4 1712 777,2 



(4.5) 



Actu ally, as a check on the algebraic manipula tion programmes (done with mathematica) that we wrote to solve 
Eq. (4.3b) we have explicitly checked that Eq. (4.5) is the unique (translation-and-rotation-invariant) solution of the 
2PN matching. 

At 3PN, we write (by using translation and rotation invariance) Y^ PN in terms of some scalar functions (here 
n Q h = (x a - x & )/r ab ; r afc = |x a - x fc |) 



Yf N (x fc , p fc ) = M a n ab + J2 N ab Pb 



(4.6) 



By imposing that the map reduces to the identity in the free-motion limit (G — * 0), it is enough to look for M a and 
N a b of the symbolic form: 

a V 2 1 In 7*12 
M ocp + — — h -5- + 

r i2 rf 2 rf 2 

ria 

where 'p™' denotes all the scalars made with pi, P2 and ni2 with homogeneity p n , i.e. 

c„ in2 „3„ 4 „ 5 (plr (PIT 2 (Pi • P2)" 3 (ni2 • Pi)" 4 (ni2 • P2)" 5 

with 2ni + 2«2 + 2n^ + 71,4 + 715 = n. We find that Yf PN a priori contains 52 unknown coefficients c„ (28 in Mi, 
12 in iVii, and 12 in Ni 2 ). We did not impose any a priori constraints on the mass dependence of the coefficients 
On (mi , 771 2) ent ering Y^ PN . (As a consequence we cannot make use of the 1 <-> 2 relabeling symmetry). Writing in 
full Eqs. ( 4.3b ) gives a linear system of 512 equations for the 2 x 52 = 104 unknown coefficients c n . In spite of this 
very high redundancy, we found that this system is compatible if and only if the arbitrary parameters tu s and A are 
related by 



A 



1987 



11 3080 

Then the solution is unique and reads (for a = 1; the solution for a — 2 being obtained by relabeling 1 <-> 2) 



(4.7) 
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Y? PN (x a , Po ) =Gm 2 < 
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where 



i (pi-p 2 )pj i (p1) 2 



m\m\ 



4 
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ni2 • Pi)(ni2 ■ gg)p| 3 (ni2 ■ P2) 2 (Pi • P2) _ 3 (nig ■ P2) 2 pj 
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mim 2 
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1 (nig ■ Pi)(n 12 ■ p 2 ) 3 1 (n 12 ■ p 2 ) 4 



r i = 167 p 2 105 ( Pl ■ p 2 ) 

1 48 m^ 16 m\rri2 
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48 mi 



(n i2 • p 2 ) 



65 TO2^ p| 25 (ni2 • pi) ; 
48 mi 
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9 (nig ■ Pi)(ni2 ■ gg) 
8 m\m2 



(4.9a) 



(4.9b) 



M 2 = 
M°1 = 
M\ = 



28387 /49 21 2 \ m 2 22 r i2 
— 7r h — In —r , 



2520 



36 32 



mi 



1 (ni2 • Pa)pf _,_ (n i2 • P2)p 2 _ 5 (n 12 ■ p 2 ) 5 
4 m\ni2 

73 (n i2 • Pi) 



24 mi 



m-2 12 m 2 

9 3 m 2 \ (ni2 • P2) 



8 2 mi 



m 2 



Y v — — - 
r i2 — 



1 (ni2 • Pi)p; 



m\m\ 



1 (ni2 
4 



P2XP1 • P2) 

mim| 



1 (ni2 

8 7 



p 2 )pi 3 (n i2 • Pi)(ni2 • p 2 ) 2 5 (ni2-p 2 ) 3 



mim| 



12 



Y 1 - 

1 12 ~ 



55 (ni2 ■ pi) 
16 



mi 
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(4.9c) 
(4.9d) 
(4.9e) 
(4.9f) 
(4.9g) 



The results ( 4.8 )— ((4.9|) give the explicit expression of the transformation (xj, p&) — > y a = Y a (xb, pj). To complete the 
knowledge of the transformation between the phase-space variables of the two descriptions one also needs the explicit 
expre ssion of the transformation (xt,,pt,) — > v a — V a (x6 , p&). This is straightforwardly obtained by inserting in Eq. 
( 4.3a ) the Hamiltonian of [0 and the results (4.4)-(fi~9|) for Y a (x&, pf,). As the explicit result is very lengthy we do 
not display it here. [Because of the availability of algebraic manipulation programmes, it is safer for the interested 
reader to rederive it directly] 

Let us mention that, as a further check, we have also tried to map the Hamiltonian H^p^(to s , LUk) containing both lo s 
and u>k to A^ armomc and that we found again that the mapping is possible only if cok = 41/24, in complete agreement 
with [ii"] . 

It should be noted that the further ambiguity parameters r[ and r' 2 present in the harmonic equations of motion 
enter our result only through some logarithmic terms in Mi [for ln(ri2/r 1 )] and M2 [for In(ri2/V 2 )]. This decoupling 
between the two particle labels (r[ entering only Yf PN , and r' 2 only Y 2 PN ) suggests (in confirmation of the discussion 
we gave above) that it is not a necessity, in the harmonic approach, to introduce the ambiguities r[ and r 2 . Indeed, 
we see that they are locally (i.e. separately for each particle) introduced by the transformation of variables between 
our (more ambiguity-free ADM result) and the variables defined by the set of prescriptions of Refs. |T(],[lj],[l3|] . We 
note also that our result confirms the finding of jl3| that r[ and r' 2 can be gauged away (by a harmonicity-preserving 
coordinate transformation). 

Anyway the most important result is that we have shown the physi cal equivalence (for invariant consequences of 
the dynamics) between the 3PN results of |^-[ll|l of those of |p^ , ^Jl7|]l^| . The invariants of the 3PN dynamics depend 
only on one ambiguity parameter, denoted tv s in the ADM work, and A in the harmonic-coordinate one. The change 
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of notation between ui s and A is given in Eq. (4.7) which agrees with the conclusion of which was restricted to 
the circular motion case. 

Could have it been different? In view of the high redundancy of the linear system we had to solve, it may seem 
quasi miraculous that the two independent results can be made to match. The compatibility we found is clearly a 
very useful check on the algebraic computations done by both groups. However, we want to point out that, as both 
groups had already checked the global Poincare invariance of their results (see for the ADM case, and [ |17|]l3|l for 
the harmonic one), the possible remaining discrepancies between the two dynamics were not very numerous. In fact 
we can count precisely the number of irreducible new coefficients entering all 3PN invariants of a Poincare-invariant 
dynamics. The simplest way to do that is to use the results of |Q on the 3PN "effective one body" dynamics JL9| . 
Using Poincare-invariance we can reduce the dynamics to that of the relative motion. Following the results of [jTsfl the 
number of irreducible new coefficients entering the relative dynamics at the rtPN level is obtained by quoticnting the 
arbitrariness in the (relative) Hamiltonian by that in a generic (relative) canonical transformation. This leaves only 
[computing the difference between Eq. (3.6) and Eq. (3.8) of [ff8|] 



(n + l)(ra + 2) 



n(n + 1) 



(4.10) 



irreducible coefficients at nPN. Moreover , one of t he c oefficients is trivial as it is given by the nPN-level expansion 
of the free-motion Hamiltonian Hq = \J m\ + + \J m\ + p\. Finally, this leaves only n non trivial irreducible 
coefficients at the nPN level, i.e., in particular, only three coefficients at 3PN. 

For instance, in the effective one body approach, these three coefficients are: 04 [the coefficient of (GM/R) A in 

in -9§M 



-Sot] 



c?3 [the coefficient of (GM/R) 3 in — ffooffiffiL an d £3 [the coefficient of P^(GM/R) 2 in the squared effective 
Hamiltonian]. These coefficients have been determined in jll|], from iJ^J^ 1 , and it was found that and 2:3 are 
unambiguously determined (independently of the oj s ambiguity) to be 



d 3 = 2(3^ - 26)v, z 3 = 2(4-3u)v, 



(4.11) 



where v = 17111712/ (mi + 1712)' 2 is the symmetric mass ratio, while 04 turns out to depend on ui s : 



a 4 



94 41 

y ~ 32' 



7T + 2w, v 



(4.12) 



In view of the sensitivity of 04 to w s , a real difference between the ADM and the harmonic dynamics could then 
only have ari sen a s possible discrepancies in the values of the only two unambiguous 3PN irreducible coefficients d% 
and 23, Eq. ( 4.11 ). As a further check, we have in fact allowed for such differences in d% and Z3 by looking for the 
matching of the harmonic equations of motion to a modified H^p^, containing two extra terms corresponding to 
variations in both c?3 and Z3. The result of this generalized matching was that the variations in da and Z3 had both 
to vanish for the matching to be possible. 



V. CONSERVED QUANTITIES AND GENERALIZED LAGRANGIAN IN HARMONIC COORDINATES 

Having explicitly obtained the transformation from ADM variables to harmonic ones which maps the two dynamics, 
we can use this map to transfer all the useful known results of the ADM approach to the harmonic one. For instance, 
Ref. has explicitly computed the ten conserved quantities of the binary system associated to global Poincare 
invariance: total energy iJ(x a ,p n ), total momentum Pi(x a ,p a ), total angular momentum Ji(x ,p a ), and the center- 
of-mass constant (boost vector) ,P° = K l (x a ,p a ) = G*(x a ,p a ) — <P l (x a ,p a ). Actually, to be able to express these 
conserved quantities within the harmonic framework one needs the inverse of the transformation (x,p) —> (y,v), i.e. 
we need to know explicitly the functions 

x Q = X a (y 6 , vb), (5.1a) 

Pa = Pa(yb,V 6 ). (5.1b) 

It is just a m atter of (somewhat involved) algebraic manipu latio ns to invert the PN-expanded map (x,p) — > (y, v) to 



get Eqs. (5.1). By straightforward insertion of the formulas (5.1), we can then (if they are needed) explicitly compute 



the following quantities in harmonic coordinates: 
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E(y b ,v b ) = H ADM (X a (y b ,v b ),-p a (y b ,v b )), (5.2a) 
Pi(yb,v b )=Y,P ai {y b ,v b ), (5.2b) 

a 

JiiYU Vb) = ^2 SikiXl{y hl v b )Pai(yb, Vb), (5.2c) 

a 

Gi{y b ,v b ) = G ADM (X a (y ,v b ),P a (y b ,v )), (5.2d) 

where the explicit expression of G ADM (x a ,p a ) is given in Ref. pj| . (Note that Gi is not conserved but satisfies 
dGi/dt = Pi. It is a useful, an d conserved, quantity in the center-of-mass frame where Pi = 0.) We shall not give here 
the explicit expressions of ( |5.2| ) as modern computer means make it safer for interested readers to do the manipulations 
themselves starting from the formulas we give. We have checked that Eq. ( |5.2a| ) agrees with the 3PN conserved energy 
obtained in Q. 

Our results allow us also to prove (without guess work) that the harmonic equations of motion derive from a 
generalized Lagrangian i harmonlc (y a , y a , Ya)) depending on positions, velocities, and accelerations. Moreover, we can 
rather simply compute the Lagrangian £ harmonlc by starting from the phase-space Lagrangian of our ADM Hamiltonian 
framework, namely 

L ADM (x a , x a , p a ) = ]T Pa • x a - H ADM (x b , p 6 ). 

a 

Indeed, it suffices (as one easily checks) to insert the transformation (x ,p ) —* (y a ,ya) 

in Z^^^ (x a , x a , p a ). This 

gives 

L harmonic (y a , y«, y„) = ]T P a (y , Yb) ■ ±a(y b , y b , y b ) - H ADM (X a (y b , y b ), P a (y b , y 6 )). (5.3) 

a 

Note that the meaning of the time derivative X l a in Eq. ([T^) is 

Therefore our constructive procedure for computing the harmonic Lagrangian automatically yields a Lagrangian 
which is linear in the accelerations y . (It was shown in |1 that it was always possible, for perturbatively expanded 
generalized Lagrangians, to reduce their acceleration dependence to be linear.) 

It should also be noted that such a linear-in-acceleration generalized Lagrangian is not unique, but is defined only 
modulo the addition of dF(y a , y a )/dt, where F(y a ,y a ) is an arbitrary^ scalar function of positions and velocities 
(only). 

When the 2PN- level generalized harmonic Lagrangian was first computed |2(],|2| , use was made of the addition of 
some i 7 2PN(ya,ya) to simplify (in a somewhat arbitrary way) the expression of £2PN(ya, Ya, Ya)- As we are not here 
playing with the addition of F we should not expect our constructive procedure ( ^.3| ) to yield a result which coincides 
with that of [p0|]2]|. We have, however, checked that our 2PN-level result is indeed equivalent to the old one, modulo 
some dF(y a ,y a )/dt. Our explicit result for the 3PN-accurate generalized harmonic Lagrangian reads (here v a = y a 
and a a = y a ) 



I , 1 . I 



L harmomc (y Q , v Q , a a ) = L N (y a , v a ) + — L 1PN (y a , v a ) + — L 2 PN(y a , v a , a a ) + — £ 3 PN(y a , v a , a a ). (5.4) 



The Newtonian and 1PN contributions to the Lagrangian (5.4) do not depend on accelerations. They equal 



10 It is, however, convenient to restrict the arbitrariness in F so as to respect the symmetries of the problem: translations, 
rotations, space parity and time reversal. 
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^ 1 2 Gmim 2 
Lvi{y a , v a ) = 2^ ^WaV a + 5 

a Z r 12 
1 22 GmiTO 2 

iiPN(ya,v a ) = -raiK) H g 

8 r l2 



3 2 7. . 1, u h . \Gm\ 

2 V 1 ~ jOl ' V 2) - 4( n 12 • v l)( n 12 ' v 2) - 



The 2PN acceleration-dependent Lagrangian L 2P n reads 

^2PN(y Q , v Q , a Q ) = L^ PN (y Q , v ) + £ 2PN (y a , v n , a Q ) + (1 <-> 2), 

where 



12 



+ (1 «-» 2). 



-^2PN 



1 2x3 , Gm 1 m 2 

= ^ m i( v i) + — = — 



'12 



Gmi 



01 | ^1 ,-02 
2PN h -°2PN 

r 12 



/ Gmi ' 
V r? 2 , 



r03 
-°2PN 



(5.5) 
(5.6) 

(5.7) 

(5.8a) 



7-01 

-^2PN 



7 27 
8 V i; 8 1 



r02 
-^2PN 



r 03 
-^2PN 



9 15 13 5 

>i • v 2 ) + ^ v ? v 2 + y ( v i ' v 2) 2 + y (n'i 2 • v!) 2 ( Vl ■ v a ) - -(nj 2 • vi)(n£ a • v 2 )v! 
-|(nj 2 • Vl )(nJ 2 • v 2 )(vi • v 2 ) - ^(n* 2 • v 2 ) 2 v 2 + ^(n h 12 ■ Vl ) 3 (nJ 2 • v 2 ) - A( n J 2 . Vl )2 (n h 2 . V2 )2 

= I v i _ I^ Vl ' V2 ^ + I v 2 + 2 (n ^ 2 ' Vl)2 ~ 2~ (n i 2 ' Vl)(n i 2 ' V2) + V^ 2 ' V2)2 ' 

1 19 m 2 

2 + Irn^' 



(5.8b) 
(5.8c) 
(5.8d) 



L 2PN = Gmim 2 ~(n 12 • v 2 )(vi • a x ) - -(n? 2 • vi)(v 2 • a x ) - -(v x • v 2 )(n5 1 2 • ai) 

+ ^v 2 ( n h 2 • ai) - i(n£ 2 ■ vi){n\ 2 • v 2 )(n^ 2 
Finally, the 3PN contribution to the Lagrangian ( |5.4| ) reads 



where 



i3PN(y a ,v a ,a a ) = L° 3PN (y a ,v a ) + i^ PN (y Q , v a , a a ) + (!<-> 2), 



Q 5 2x4 , Gm 1 m 2 

L3PN= 128 mi(Vl) + ^2~ 



rOl 

-^PN 



~ h - -^PN 
'12 



' 12 



V" 



'12 



1-04 
-^3PN 



-^3PN = Gmim 2 I L 3PN H — L 3PN 



rOl 
-^3PN 



'12 



^(V?) 3 - g (V 2 ) 2 ( V1 • V 2 ) + |(V 2 ) 2 V 2 + f V 2 ( V1 • 



(5.8e) 

(5.9) 

(5.10a) 
(5.10b) 



\ 2 2 / n2 *^ / x" 3 

v 2) - ggVfCVi -V 2 )V 2 - — (Vi • V 2 ) 



-g(n^ 2 • vi)(4, • v 2 )(v?) 2 - |(n? 2 ■ v 2 ) 2 (v 2 ) 2 + g(nf 2 ■ Vl ) 2 v 2 ( Vl • v 2 ) 

1,1, \/h \ 2 / \ 11 , \22/ \ *^ / h \ 2 2 2 

-^( n i2 ' v i)( n i2 ' v 2 )v 1 (vi • v 2 ) + y(ni 2 ■ V 2 ) Vi(Vl • v 2 ) - ^( n i2 " V l) V 1 V 2 
+^( n i2 • v i)( n i2 • v 2 )v 2 v 2 - ~(nf a ■ Vl ) 2 ( Vl • v 2 ) 2 - ^(n' 1 , • vi)(nf a • v 2 )( Vl • v 2 ) 2 
+ ^K 2 • V!) 3 (nh 2 • v 2 )v? + A (n h 2 . Vl) 2 (n h 2 . V2) 2 V 2 _ ^ (n h 2 . Vi)(n h 2 . V2) 3 V 2 + 19 (n h 2 . V2) 4 V 2 
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~(*12 ■ vi) 4 (vi • v 2 ) - (n h 12 • Vl ) 3 (n? 2 • V2 )(vi • v 2 ) + |(n^ 2 • Vl ) 2 (nJ 2 • v 2 ) 2 ( Vl • v 2 ) 
- A(n? 2 • vi) B (n^ 2 • v 2 ) - A (n h 2 . Vl) 4 (n h 2 . Va)2 + 15 (n h 2 . Vi) 3 (n h 2 . ^ 

Y (V 2 ) 2 ~ ^ V '( V l ' V 2) + f V ' V 2 + • V 2) 2 - ^(Vl • V 2 )V 2 + ^(v 2 ) 2 - ^(nj 2 • Vl ) 2 V 2 



(5.10c) 



3 
44 



+19(n? 2 • Vl )(n? 2 • v 2 )v 2 - ^(n 1 ^ • v 2 ) 2 v 2 - \(n\ 2 ■ Vl ) 2 ( Vl ■ v 2 ) - ^(n h 12 ■ Vl )(n? 2 • v 2 )(vi • v 2 ) 
+ ^(n? 2 • v 2 ) 2 ( Vl • v 2 ) + ^(n* 2 • Vl ) 2 v 2 - H(„h a . Vl )(n? 2 ■ v 2 )v 2 + l(n? 2 • v 2 ) 2 v 2 - H (n?2 . Vl )4 
+f(n h 12 ■ Vl ) 3 (n^ 2 • v 2 ) - ^(n* 2 • vi) 2 ^ • v 2 ) 2 + ^ (n 1 ^ • Vl )(n? 2 • v 2 ) 3 , 



15611 + I 41 ^ 2 305^ m 2 



1260 + V 128 



144 



mi 



17501 , (AI _ 2 439 



1260 V 64 144 



/ll 2 io» \ m 2 



mi 



5 / 41 2 305\ m 2 l 2 T8243 /123 2 383 \ m 2 

7T V, — h 7T 

128 144/ mi J 2 [210 \ 128 48 / m i 



(vi • v 2 ) 
(n^-vi) 2 



15541 /123 2 889 \ m 2 
f ~48~ y mrj~ 



420 
22 



(ni2-vi)(nj 2 -v 2 ) + 



383 123 2 \ m 2 



48 128 J mi 



(n? 2 -v 2 ) 2 



22 22 

yv? + y(vi • v 2 ) + 22(n'i 1 2 • vi) 2 - 22(ni 1 2 • vi)(nj 2 • v 2 ) 



3 /9707 \ m 2 22 m 2 



V 420 



- u> 



H0 2 W UU 2 I 1 2 

1- -7. In—, 

mi 3 mi r\ 



(5.10d) 



(5.10e) 
(5.10f) 



15 5 7 

- — ( n 'i2 ' v i)( v i • v 2)(vi • ai) + -(n 12 • vi)v|(vi • a x ) + -(n 12 • v 2 )v 2 (vi • ai) 

1 5 5 

-2( n 'i2 ' v 2)( v i ' v 2)(vi • ai) - -(n^ 2 • v 2 )v 2 (vi • ai) - -{n\ 2 ■ vi) 2 ^ • v 2 )(vi • ai) 

3 5 15 

-^( n i2 ' v i)( n i2 ' v 2 ) 2 (vi • ai) + ^(n? 2 ' v 2 ) 3 (vi • ai) - — (vl\ 2 ■ vi)v 2 (v 2 • ai) 

+ ^( n i2 ' vi)(vi • v 2 )(v 2 • ai) - j(n^ 2 • v 2 )v 2 (v 2 • a x ) - ^(n^ • v 2 )(vi • v 2 )(v 2 • a x ) 

3 5 3 

-^( n i2 ■ vi)v 2 (v 2 • ai) + ^(n^ ■ vi) 3 (v 2 • ai) + -(n 1 ^ • vi) 2 (n^ 2 • v 2 )(v 2 • ai) 

3 15 5 

-g( n i2 ' v i)( n i2 ' v 2 ) 2 (v 2 • ai) - — v 2 (vi • v 2 )(n5 1 2 • ai) + -vfv^n^ • ai) 

+l(vi • v 2 f(n\ 2 ■ ai) - ^(vi • v 2 )v 2 (n 1 1 1 2 • ai ) + ^(v 2 ) 2 ^ • &1 ) 

-\{tl\ 2 • vi)K 2 • v 2 )v 2 K 2 • ai) - |(n? 2 ■ v 2 ) 2 v 2 (n^ 2 • ai ) + \(n\ 2 ■ vi) 2 (vi • v 2 )(n? 2 • &1 ) 

3 3 5 

+ 2( n 'i2 ' v i)( n i2 ' v 2)(vi • v 2 )(n 12 • ai) - -(n 12 • v 2 ) 2 (vi • v 2 )(n 12 • ai) - -(n 1 ^ • Vi) 2 v§(n 12 ■ ai) 
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+~(n* 2 • Vl )(n? a • v 3 )vl(nj a • a x ) + ~(n> 2 • Vl ) 3 (n? 2 • v 2 )(n* 2 • a x ) + ~(n* 2 ■ Vl ) 2 (n* 2 • v 2 ) 2 (n^ • a x ) 
-^(nJa-ViJCnJa.vaJVJa-ai), (5.10g) 

12 19, h \/ \ 185 , h \/ \ 205 , h \/ \ 2/ h \ 

L 3PN = ~^j;( n 12 • V l)( V l • & l) ~ ^("12 ' V2 )( Vl ' &1 ) + "^"^ ' Vl )( V 2 ' &1 ) + y V l( n 12 ' a l ) 

~^2"(vi ' vs)(nf 2 • ai) + ^v 2 (n^ 2 • ai) + ^(n^ 2 • vi) 2 ^ • ai) - y(n? 2 • vi)(n£ 2 • v 2 )(n^ 2 • ai) 

5.i, . , . 235 /V. \ * s. 31 . u , . . 35 9 / h \ 
-|( n i2 ' v i)( v i ' a 2j + -24"( n i2 • v 2)(vi • a 2 ) - y(nf 2 • Vi)(v 2 • a 2 ) - jvf(< 2 ■ a 2 ) 

+ff(vi • v 2 )(n? 2 ■ a 2 ) - ^« 2 ■ Vl ) 2 (n* 2 • a 2 ) + y (n 1 ^ ■ Vl )(n* 2 • v 2 )(n* 2 • a 2 ). (5.10h) 



VI. CONCLUSIONS 

To conclude, let us emphasize that the equivalence, established in this paper, between the existing independent 
approaches to 3PN dynamics is important because it confirms the basic soundness of both approaches. It shows that 
the quite different regularization procedures devised by the two groups are physically equivalent. None can claim 
to be mathematically 'better' or 'more correct' than the other one. We have, however, pointed out that the ADM 
regularization: (i) is significantly simpler to define and apply in practice, (ii) leads to the introduction of a minimal 
set of regularization ambiguities (without extra gauge-related ambiguities). 

So much for the good news. The bad news is that having proven the physical equivalence between the two approaches 
sheds no light on the problem of the 'static ambiguity' u) a . In fact, it is sobering to note that the enormous work 



which went into both 3PN investigations § M^JlJj|^Jrj,y^y7|,|3| , and which led to the explicit evaluation of 0(100) 
3PN coefficients (not to mention the 0(lO^intermediate expressions which had to be computed and manipulated) 
succeeded in getting only two out of the three irreducible 3PN coefficients mentioned above (^3 and Z3; with 04 staying 
ambiguous). This is all the more a pity that it was shown in |l8| that the 3PN- level predictions for the physically 
most important quantities (dynamical behaviour [ fl8[ , and gravitational wave emission pl|] near the transition between 
inspiral and plunge of a binary black hole) vary quite significantly when tu s is allowed to vary within the plausible 
range of —10 < uj s < 10. This makes it urgent to further clarify the origin of the 'static ambiguity'. 

Roughly speaking, it seems that this ambiguity is due to the breakdown, at 3PN, of the possibility of modelling 
extended (compact) objects (neutron stars or black holes) by delta-function sources. This is somewhat surprising 
because it has been shown long ago that the extended nature of the objects (violation of the 'effacing principle') 
should show up only at 5PN (see Ref. p. 86). One must probably use new techniques (or extend to 3PN-level 
old techniques, such as the matchin g te chnique used in ||) to solve this problem.^ We note that it would be nice 
if the 'effective one body' approach [|19|, which is so efficient in condensing the invariant content of the dynamics to 
a few coefficients, could be developed into a calculational technique, giving explicit algorithmic recipes for directly 
computing the three 3PN irreducible coefficients. 
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Let us recall that Ref. jL8| found several independent arguments suggesting that ui s ~ —9, and maybe ui s = — + |j7r 2 . 
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